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Abstract 

It is shown that computer simulations can qualitatively reproduce ex- 
periments, where a powder of cohesive, round, hard particles is periodically 
deformed at constant volume. Two types of initial configurations are consid- 
ered: Uniaxially precompacted ballistic deposits and biaxially precompacted 
DLA-clusters. Both initial configurations had the same volume fraction, but 
due to the different precompaction procedure completely different principal 
stresses. After a transient which lasts only less than a period, the stresses 
follow the same periodic function, i.e. the powder forgot its history. 

1. INTRODUCTION 

One of the fundamental open questions in the physics of powders is how to 
characterize their state uniquely, the particle configuration in general depending 
on the deformation history of the powder. This is in marked contrast to ther- 
mal equilibrium systems, which statistically explore all available configurations, 
so that temporal correlations normally die out quickly. A powder, however, is 
"jammed" or frozen in, if it is not driven by external forces. Even if it flows, 
it does not explore its alternative configurations on a time scale which is short 
compared to the inverse strain rate. 

This has macroscopic consequences. For example the stress in a compacted 
cohesive powder not only depends on the solid fraction but also, by what kind of 
compaction the solid fraction was reached. It is believed, that at least the fabric 
tensor is needed in order to determine the stress state. Any microscopic theory of 
the dynamical evolution of the state and hence of powder flow requires an answer 
to the question, how to characterize such a state. Therefore it is very important 
also for macroscopic predictions needed e.g. in process engineering. 
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In this paper we address only a partial aspect of this question. It is known, that 
the stress components in a powder, which is quasi-statically deformed at fixed vol- 
ume, approach values independent of the precompaction history lH). This stress 
state depends on the volume fraction, but also on the total deformation at constant 
volume: As one continues to shear, the stress components increase. The point 
we want to make in this paper is that there are also other ways to make the pow- 
der forget its precompaction history. As an example we investigate the periodic 
kneading of the powder at constant volume. In contrast to steady state flow the 
principal axes switch back and forth periodically in this case. This has the advan- 
tage that the sample returns to its initial shape periodically. In every period steady 
state flow properties may be reached or not, depending on how long the period 
lasts. 

Periodic kneading was studied in a different context previously Q. There it 
was shown for limestone that the stresses are periodic functions with decreasing 
amplitude. After a few periods the amplitude becomes constant. In the present pa- 
per we compare experiments on carbonyl-iron powder and computer simulations. 
In contrast to limestone the particles are spherical and very hard. As we are going 
to show, the amplitude of the periodic stresses becomes constant already within 
one cycle in this case. Simultaneously all stress tensor differences due to different 
sample preparations vanish within this short transient, i.e., the powder forgets its 
history. 

2. METHODS 

2.1 Experimental Setup 

The true biaxial shear tester allows the deformation of a brick shaped bulk 
solid specimen. The deformation mechanism is composed of a bottom and a top 
plate which are fixed at a distance of 36.5 mm from each other and four side 
plates which can be moved independently from a maximum distance of 130 mm 
to a minimum distance of 70 mm, always staying perpendicular to their adjacent 
ones. Thus the bulk solid specimen remains brick shaped and principal strains 
equal the normal strains (Fig. [it): 

eij = 0, where i j; eu = ei . (1) 

The side plate's velocities can be set in steps of 0.03 fimJs up to a maximum 
velocity of 65 /im/s corresponding to maximum shear rates of 0.1% per second. 
On the specimen's border, the complete stress state is measured with three five 
component load cells |l3|, which are installed in the bottom and in two perpendic- 
ular load plates, LC 1, 2, and 3 in Fig.[lj5. For the tests with the true biaxial shear 
tester, silicon grease is spread on bottom, lateral and top plates. Then they are cov- 
ered with highly flexible rubber membranes in such a way that the membranes are 
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Figure 1: a) Schematic view of the sample chamber explaining the axes, b) The 
true biaxial shear tester. 

evenly deformed with the borders of the powder specimen. The friction between 
specimen and load plates thus is minimized, hence shear stresses are negligible. 
Therefore, the normal stresses on the specimen borders are principal stresses. The 
principal stress axes coincide with the deformation axes (Fig.[lt): 

aij = 0, where i j; an = ai. (2) 

In this paper we use the engineering sign convention for the stress, i.e. a com- 
pressive stress is defined as positive. Deformations and stresses are continuously 
measured using a process computer, which also controls the deformation process 
lis. As experiments with finite strains are carried out, we define the strains using 
the natural strain increment t4J and obtain 

6. = ln^ (3) 

for the strains in x-, y- and z-direction and for the volume strain respectively: 

ey = In — = In - — - — — = + ey + (4) 

In order to perform tests at constant volumes, one can easily control the condition 
+ =const. (e^ = during the complete measurement due to the fixed bottom 
and top plates). For test preparation, the powder is first sieved in z-direction into 
the tester. Then, it is precompacted in the same direction in order to homogenize 
the sample. After this, the sample is sheared with = —ty, i.e., the sample's 
volume remaining constant. The directions of shear are changed periodically. All 
experimental results presented in this paper have been obtained with carbonyl-iron 
powder with a median particle size of about 2 microns. This powder was chosen 
because its primary particles are perfectly round, in analogy to the particles used 
in the computer simulations. 
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2.2 Simulation Method 

For the computer simulations we use a discrete element method, the Contact 
dynamics method 0. In contrast to soft particle molecular dynamics, the volume 
exclusion of perfectly rigid particles and Coulomb's friction law are implemented 
exactly [6|. The cohesion between the particles is included in a simplified model 
as an attractive force acting within a range leading to a force and an energy barrier 
to break a contact between two particles |El[8l|9|. Furthermore, rolling friction 
between the round particles is implemented, also leading to the stabilization of 
large pores EIHl. We use two dimensional simulation because the influence of 
the torsional degree of freedom for a contact present in three dimensions is not 
well understood up to now. 

All simulation results are given in the following units: As unit of length we 
take the particle radius r, as unit of mass pr^, where p is the mass density, and as 
unit of time the inverse strain rate je^l"^ = \^y\^^ (which was constant during the 
periodic kneading). In these units the simulation parameters were: The cohesion 
force between two particles, = lO^^pr'^ l^xP^ the range of the cohesion force, 
d = 10^^ r, the friction coefficient, p = 0.3, and the rolling friction coefficient, 
/iR = 0.2 r (which is the torque resisting the rolling motion divided by the normal 
force Fn + Fc at the contact). We like to point out that for the typical values of the 
carbonyl-iron powder experiment [r ^ 1 /im, p ^ 10^ kg/m^ and le^;] ~ 10~^ s~^] 
our simulation parameter Fc ~ 10~^ N is a reasonable value for the cohesion force 
between the particles. 

The walls were cohesive, but frictionless, so that no tangential stresses could 
occur, as in the experiment. This guarantees, that the principal axes of the stress 
and the strain tensors were coUinear. 

3. RESULTS 

For the two-dimensional simulations, we prepared two different kinds of initial 
configurations in the following way (cf. Fig. |2l): One is a ballistic deposit ifTOl . 
which was precompacted in the direction of the deposition [8l up to a volume 
fraction of u = 0.42. The other is an off-lattice DLA cluster [llj centered in a 
square box, and then precompacted biaxially to the same volume fraction, shrink- 
ing the square box symmetrically. Note that the first preparation procedure tends 
to lead to anisotropic, but homogeneous configurations, while the second one re- 
sults in inhomogeneous, but isotropic configurations (the initial DLA cluster is 
essentially round and has a higher density at the center). This was done delib- 
erately to demonstrate that even in this case both configurations reach the same 
(periodic) stress/strain relationship. 

The precompacted systems are now used as initial configurations for the pe- 
riodic shearing at constant volume. As in the experiment, the shearing was per- 
formed with constant logarithmic strain rate. The maximal strain was 36% (with 
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Figure 2: a) Initial configuration 1: Ballistic deposit, b) Initial configuration 2: 
DLA cluster, c) Configuration 1 precompacted uniaxially in the deposition direc- 
tion, d) Configuration 2 precompacted biaxially (keeping the square symmetry), 
e) Configuration 1 after one cycle, f) Configuration 2 after one cycle. 

respect to the square shape). Figure |3] shows the time dependence of the two 
principal stresses, averaged over 20 independent runs (with different, but simi- 
larly prepared initial configurations), including the precompaction. Moreover the 
data were smoothed by averaging them over a running time window of 5% of the 
period. 

Clearly, the stresses are qualitatively different during the precompaction: For 
configuration type 1 (ballistic deposit) and uniaxial precompaction one gets aniso- 
tropic stresses, the larger stress belonging to the direction of compaction. For con- 
figuration type 2 (DLA-cluster) and biaxial precompaction both principal stresses 
increase in the same way, but when the square is deformed into a rectangular shape 
at constant volume (first 1/4 cycle), the principal stress in the expanding direction 
drops to a small positive value, while the other component stays large. This be- 
havior agrees qualitatively with the experiment with the biaxial precompaction. 

Fig. a 

The time dependent stresses within the first half cycle (t/tcycio = to 1/2) in 
Fig.lSK are qualitatively different from the ones in all subsequent half cycles. Sim- 
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Figure 3: Simulation results: Principal stresses as functions of time during peri- 
odic deformation a) for configuration type 1 [volume constant for t/tcyde > 0], 
and b) for configuration type 2 [volume constant for t/tcyde > —1/4]. 



25 



■^01 5 
CL 



1 


1 1 1 


1 1 1 1 1 1 1 


1 1 1 

.. a 


1 








X 

a 












y 

- a 














z 








d 


I — 




; ' / / 




1 




■■■■>■■■ 




, J I.- 


— \ — (....I.-- 


— 1 



5000 10000 15000 

t[s] 



Figure 4: Experimental results: Principal stresses as functions of time during 
periodic deformation of biaxially precompacted carbonyl-iron powder. 
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ilarly the first 1/4 cycle (t/tcycic = —1/4 to 0) in Fig. |3j5 is different. However, 
after these transients, it seems that the same stress pattern is repeated period after 
period, also in the experiment Fig. SI In order to compare the stress patterns for 
the two cases. Fig. [3^ and b, we averaged them over 9 repetitions. The results are 
shown in Fig|5l Within the error bars the stress pattern for the samples obtained 
from uniaxially precompacted ballistic deposits is the same as for those obtained 
from biaxially precompacted DLA-clusters. Moreover Oy simply coincides with 
ax shifted by half a period, as it should. This shows, that after the transients the 
stresses do not show any dependence on the preparation history. This is also re- 
flected in Fig|2t and f, which show that after one cycle the particle configurations 
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Figure 5: One period of cr^. and of Oy, averaged over all repetitions after the initial 
transient, a) For v = 0.42 as in Fig.|3l b) for u = 0.63. In both cases the curves for 
initial configurations of type 1 [full line] respectively type 2 [dashed line] coincide 
within the statistical errors. 



1 1 1 1 1 1 1 1 

— ballistic (uniaxial) 
DLA (biaxial) 








^ I.I, 




7 



are qualitatively indistinguishable. 

Figure |5j5 shows the stress pattern for a higher solid fraction. Here the max- 
imum strain compared to the square shape was 16%. Again one gets the same 
pattern in spite of the different types of original configurations and precompaction 
histories. However, the shape and amplitude of this pattern is completely different. 
This shows that the periodic stress pattern depends on the solid fraction. 

4. CONCLUSIONS 

We investigated the strain-dependent stress in a cohesive powder, which was 
periodically deformed at constant volume and strain rate. It was shown that con- 
tact dynamics simulations agree qualitatively with experiments on carbonyl-iron 
powder. 

For the simulations two types of initial configurations were prepared. One 
had anisotropic stress and fabric tensors (uniaxially compacted ballistic deposits), 
the other was isotropic (biaxially compacted DLA-clusters), but both types of 
configurations had the same volume fraction. These initial configurations were 
then periodically deformed at constant volume. After a transient which lasts only 
less than a period, one obtains for both types of initial configurations the same 
periodic stress pattern. This leads us to the conclusion that periodic shear like 
steady state flow makes the powder forget its history quickly. It was also shown, 
that the periodic stress pattern depends on the solid fraction. 
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6. NOMENCLATURE 

e - strain tensor, dimensionless 

a - stress tensor. Pa 

L - size of the sample chamber, m 

V - volume of the sample chamber, m^ 

r - particle radius, m 

p - mass density, kg m~^ 

d - range of the cohesion force, m 

Fc - cohesion force, N 

Fn - normal force, N 

fj, - coefficient of tangential friction, dimensionless 
//R - coefficient of rolling friction, m 
u - volume fraction, dimensionless 
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